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Plots  of  acoustic  power  versus  frequency,  and  acoustic  pressure  direc¬ 
tivity  at  selected  frequencies,  are  given  for  the  separate  cases  of  thin, 
thick  and  exact  plate  theories  in  which  the  excitation  is  a  vertical  time- 
harmonic  point  force;  they  illustrate  the  range  of  applicability  of  the 
approximate  theories.  Plots  of  acoustic  power  and  pressure  directivity  are 
given  for  exact  plate  theory  in  Which  the  excitation  is  a  longitudinal  force. 
Plots  of  pressure  transmission  and  reflection  coefficients  axe  given  for  exact 
plate  theory.  Contour  and  perspective  plots  of  pressure-frequency-angle  data, 
obtained  from  exact  theory,  provide  vivid  illustrations  of  the  acoustic 
properties  of  the  plate. 
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Pestell  and  Janes  [1]  have  calculated  faur-field  sound  radiation  frost 
layered  media,  excited  by  tine-harmonic  forces,  by  the  method  of  dynamic 
stiffness  coupling  in  Which  each  element  of  the  system  is  represented  by  a 
'spectral*  dynamic  stiffness  matrix.  Based  on  this  work,  Spicer  [2]  has 
simplified  the  mathematical  derivation  of  the  stiffness  matrices  and  has 
included  uniform  mean  subsonic  flows  in  the  acoustic  fluid  elements.  Fortran 
programs  arising  from  the  aforementioned  studies  have  been  used  to  predict 
sound  radiation  from,  and  free-wave  propagation  in,  a  variety  of  layered 
systems.  Additionally,  a  Fortran  program  has  been  written  to  calculate 
reflection  and  transmission  coefficients  due  to  plane-wave  excitation. 

The  acoustics  of  fluid-loaded  infinite  plates  has  been  studied  exhaus¬ 
tively  in  the  literature,  but  only  a  restricted  number  of  graphical  presen¬ 
tations  have  been  published.  The  mathematics  of  the  'thin'  and  'thick'  plate 
theories  can  be  found  in  the  standard  text  on  fluid-structure  interaction  [3]. 
Evseev  et  al  [4]  and  Sergeev  [5]  have  studied  sound  radiation  using  plate 
theories  which  include  the  effects  of  longitudinal  vibrations  as  formulated  by 
Lyamshev.  Freedman  [6]  gives  a  limited  comparison  of  the  transmission 
coefficient  obtained  from  exact  and  approximate  plate  theories  and  discusses 
the  effect  of  including  the  first  symmetric  longitudinal  wave  motion  in  the 
thin  plate  model.  Fiorito  et  al  [7]  have  used  a  resonance  formalism  theory 
which  provides  physical  insight  into  the  peaks  and  dips  that  occur  in  plots  of 
transmission  and  reflection  coefficients  versus  angle  of  incidence.  Woolley 
[8]  has  extended  thick  plate  theory  to  include  symmetric  modes  of  vibration. 

It  is  the  purpose  of  this  report  to  extend  the  range  of  published  numeri¬ 
cal  results  of  sound  radiation,  reflection  and  transmission  of  an  Infinite 
steel  plate  which  is  iammrsed  in  water.  Subsequent  reports  will  include 
numerical  results  of  the  acoustics  of  layered  madia  and  will  also  discuss 
wavenumber  versus  frequency  plots  obtained  from  the  dispersion  relation  of  the 
layered  system.  The  maserical  results  contained  herein  will  help  to  illustrate 
and  aid  understanding  of  the  acoustics  of  a  steel  plate  t  they  will  also 
provide  reference  plots  for  comparisons  with  plots  obtained  from  simplified 
theories. 

Section  2,  below,  gives  some  details  of  the  mathematics  needed  for  both 
exact  and  approximate  plate  theories,  the  exact  theory  being  based  on  the 
multilayer  theory  of  Spicer  [2].  The  numerical  results  are  discussed  in 
Section  3  and  they  are  followed  by  some  concluding  remarks  in  Section  4. 


2.  HMTOMICS  of  pboblem 

<a)  Trtn— Udiifln  wieeei™  Bast 

In  Figure  1A  a  time  harmonic  plane-wave  is  incident  on  a  layered  system 
Whose  elements  are  governed  by  the  equations  of  acoustics,  elastodynamics  and 
visco-dynamics.  The  incident,  reflected  and  transmitted  waves  must  be  of  the 
form 
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Pt  -  Pt«xp<  iax+itfy-iy^ ) 


(2.1) 


Pr  -  PjR.exp^ax+ijSy+iy^) 

Pt  ”  PiT.exp< iax+i^y-iy2[z+H) ) 


«h«n  Pj^  is  the  amplitude  of  the  incident  wave;  R  is  defined  as  the  complex 
pressure  reflection  coefficient;  T  is  defined  as  the  complex  pressure  trans¬ 
mitted  coefficient;  a,  0,  y*  and  y2  are  wave-numbers  which  must  satisfy  the 
constraints 


yi“  +v<*i"«2-02> 
y2-  +v(*l-a2-02) 


(2.2) 


because  the  pressures  must  satisfy  the  acoustic  wave  equations  in  their  re¬ 
spective  domains;  k1  and  k2  are  acoustic  wavenumbers ,  u/cl  and  u/c2;  the 
time  harmonic  factor  exp( -iwt )  is  osd-tted  throughout. 
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Spicer  [2]  has  presented  the  mathematics  needed  to  generate  the  matrix 
relation  connecting  interface  displacements  and  forces,  viz. 


[Z(o,0)][O(a,0)]  -  [E(a,0)]  (2.3)  i 

N*N  11*1  11*1  || 

r 

where  M-3N+3  with  N  being  the  number  of  layers;  [Z(a,0)]  is  the  system 
dynamic  stiffness  matrix; 


[U(a,0)]  -  [uxl,uyl,u2l,  . 


,u  ,u  ,u  ] 
'  xn  yn  znJ 


is  a  column  vector  of  displacements  at  the  n  layer  interfaces,  where  n-N+1; 

[E(a,0)]  -  [0,0, -ZP^ . 0,0, 0]T 

is  a  column  vector  of  interface  stresses,  assumed  to  be  positive  When  acting 
in  the  positive  x-direction,  that  takes  the  above  form  for  plane -wave 
excitation  because  the  excitation  is  simply  the  ’blocked  pressure’  [3], 
-tPi+PrJz-0*  factor  exp(  iox+i0y )  has  been  omitted  from  both  sides  of  the 
equation. 

The  normal  displacements  of  the  upper  and  lower  surfaces  of  the  system, 

vis. , 


U  (a,0)  -  u  (a,0)exp(iaxti0y) 
1  zi 

0  (a,0)  -  u  ( a, 0 )exp< iax+10y ) 
z  zn 


(2.4) 
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are  found  by  matrix  inversion  of  equation  (2.3) 


'  *  .-wt ~  %  •  -«  \ 


The  boundary  condition 

[dp±/»z  +  dp t/9z}zm0  - 

that  requires  continuity  of  fluid  particle  and  interface  displacements  at  the 
upper  surface  gives 


R  -  1  +  p^  uel(a,/3)/iylPi 


(2.5) 


and  a  similar  boundary  condition  applied  at  z— H  gives 


T  -  -p2«  Uzn(af0)/iy2P. 


(2.6) 


Por  the  particular  case  shown  in  Figure  1,  the  wavenumbers  and  angle  of 
incidence,  eif  are  related  by  the  simple  equations 

a  «  k1sin(ei) 

0-0 

yt  -  k1cos(0i) 

y2  “  +V(k2"ltt8in2(0i)l 

The  angle  of  transmission ,  et,  can  be  found  from  Snell's  law,  viz. 


which  gives 


X2sin(et)  -  k1sin(ei) 


-1  2  2  2  1/2 
«t  -  tan  [Xlsin©i/(k2-X1sin  8^  ] 


(2.7) 


(2.8) 


It  may  be  deduced  from  equation  (2.8)  that  when  c1>cz,  there  is  always  a  real 
angle  of  transmission)  and  then  c1<c2  there  is  a  critical  angle  of  inci¬ 
dence,  0^-sin-1(c1/c2 ),  beyond  which  there  is  no  transmitted  wave.  This 
phenomenon  is  discussed  in  numerous  texts,  for  example  (9]. 


The  pressure  reflection  and  transmission  coefficients  in  decibels  are 
defined  as 


Rd  -  20.0  x  log10IR| 


(2.9) 


T.  “  20.0  x  log,_|T| 


WHWP  Fa  PJ  r: 


(b)  Sound  Radiation. _ Enact 


In  Figure  2B,  time-harmonic  point  forces  and  sources  excite  the  layered 
system,  whose  fax-field  pressure  in  the  upper  half-space  is 

p<R,d,*)  -  -p1u2uzl(ot,0)exp<ik1R)/2irR  (2.10) 

in  which  the  'stationery  phase*  [3]  wavenumbers  are 


a  «  k1sin(0)cos(4>) 
0  -  k18in(d)sin(0) 


(2.11) 


The  displacement  uzl(a,/3)  is  found  by  matrix  inversion  of  equation 
(2.3)  in  which 


(E(<*,0)]  - 

A  vertical  point  force,  P0, 
represented  by,  for  example 


. 'xn'V'xn*®*.1 


located  on  the  upper  interface  at  (x-xo,y-y0 ), 


is 


*zf  Vx*"laVi*r0>  (2‘12) 

and  a  point  source  of  sound,  p0axp< ik2R)/R,  located  at  distance  z0  below  the 
bottom  interface  gives 

S*n  "  4ffiP0«‘P<-i«o_i^yO+i,*Oly2>/y2  <2*13> 

The  sound  level  in  dB  ref.  l  p*t/m2  at  lm  is  defined  as 

dB  -  20.0  X  logiQ|p(R-l,e,0)|  +  120.0  (2.14) 

For  the  case  of  point  force  excitation  with  peak  amplitude  F0,  the  decibels 
as  defined  above  must  be  regarded  as  for  force  P0  rats. 


(c)  Sound  Radiation  of  Plata.  m»t-a 

For  the  special  cases  of  thin  and  thick  ( Timoshenko-Hind 1 in )  plate 
theories  the  sound  pressure  is  again  given  by  equation  (2.10)  in  which,  for 
point  force  excitation 


u8l(o,0)  -  FQexp(-laxo-10yo)/Z(a,0) 


where 


or 


Z(a,0)  »  D(o2+02)2-w2p<h-ip1«2/y1-ip2w2/y2 


(2.15) 


(2.16) 


n 


f 
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Z(a,/3)  -  {[(D<a2+/32)-psh3u2/12)(a2+^2-psM2/KG)-«2psh] 

/[l+0<  a2+/32  )/KGh-p#h2«2/12KG]  J-ip^/y^i  P^/yz 


respectively .  o  la  the  plate  flexural  rigidity;  p8  la  its  density  and  h 
its  thickness;  G  is  the  shear  modulus  and  X  *  rr2/ 12  is  a  shear  correction 
factor;  yt  and  y 2  are  given  by  equation  (2.2).  The  plate  theories  are 
discussed  in  some  detail  in  the  standard  text  [3]  on  fluid-structure 
interaction. 


The  radiated  power,  in  the  upper  half-space,  is  defined  as 


P  -  ( 1/  pc ) 


2  n  n/2 

i  I  ' 


0  0 


|p(R,e,4>)|2  R2sin( e )ded4> 


which  becomes,  on  using  equation  (2.10) 

2ir  n/2 


(pM4/4ir2c)  |u  (a, /3)|  sln( 0 )ded4> 

zi 


0  0 


(2.17) 


(2.18) 


For  the  special  case  of  axisymmetric  excitation,  it  is  relatively  easy  to  show 
that 


(pw4/ 2»rc)  |u  ,(<x->k  sine) | 2  sin(d)de 

I  Zl  1 


(2.19) 


The  sound  power  level  in  dB  ref  lpw  is  defined  as 


dB  -  10.0  x  logiQ|P|  +  120.0 


(2.20) 


For  the  case  of  point  force  excitation  with  peak  amplitude  F0,  this 
definition  of  power  level  must  be  regarded  as  for  force  F0  ms;  hence,  the 
customary  multiplication  factor  of  one-half  has  been  omitted  in  equation 
(2.17). 


In  Figures  2-17,  which  shew  transmission,  reflection  and  radiation 
levels,  the  following  SI  const  amts  were  used; 


ste«l  plate i  x 


■  10. 437x10 10  M  -  7. 558x10 10 
Z  -  19.5x10*°  <T  -  0.29 
pg  -  7000.0  h  -  0.05 

water:  p  -  iooo.o  c  -  isoo.o 

Plate  damping  was  Included  In  the  computations  by  setting  the  elastic 
constants  to  complex  values,  for  example  EsE(l-in),  Where  ti  is  called  the 
hysteretic  loss  factor. 


(b)  Acoustic  Power 

Figures  2  and  3  show  acoustic  power  radiated  into  the  upper  half-space  by 
a  plate  that  is  excited  by  a  lkN  point  force.  The  lower  half-space  contains  a 
vacuum.  The  power  levels  were  Obtained  via  numerical  integration  of  equation 
(2.19).  Small  irregularities  in  the  plots  axe  due  to  the  low,  but  adequate, 
accuracy  of  the  integrations. 

Figures  2A,  2B  and  3A  form  a  comparison  of  power  radiated  by  thin,  thick 
and  exact  plate  theories  in  which  the  excitation  is  applied  at  the  lower  sur¬ 
face  in  the  vertical  direction.  The  effects  of  increasing  the  damping  factor 
from  0.01  to  0.1  are  quantitatively  similar  on  all  three  plots,  viz.,  a  very 
small  effect  at  frequencies  below  5kHz  rising  to  an  8dB  reduction  at  50kHz. 

The  levels,  however,  depend  on  the  particular  plate  theory  used;  the  thin 
plate  theory  is  inadequate  above  5kHz  and  the  thick  plate  theory  is  wanting 
above  40kHz. 

A  plot  of  acoustic  power  radiated  by  a  longitudinal  force  applied  to  the 
lower  surface  is  given  in  Figure  38.  A  comparison  of  Figures  3A  and  3B  shows 
that  while  the  power  radiated  by  a  longitudinal  force  is  significantly  less  at 
frequencies  below  10kHz,  it  is  only  about  6dB  less  than  the  power  radiated  by 
a  vertical  force  at  higher  frequencies.  Because  this  plot  was  also  obtained 
via  equation  (2.19),  allowance  had  to  be  made  for  the  fact  that  the  pressure 
varies  in  the  azimuthal  direction  as  cos($),  i.e.  this  equation  was  multiplied 
by  one-half. 


(C)  Radiated  PX9BBUXQ ■ _ comparison  of  Theories 

Figures  4-9  show  the  variation  with  angle  9  of  the  fax-field  pressure, 
at  selected  frequencies,  of  the  three  plate  theories.  The  lower  half-space 
contains  a  vacuum.  The  excitation  is  a  lkN  point  force  which  is  applied 
vertically  at  the  lower  surface. 

At  2kHz,  Figure  4,  the  sound  levels  of  the  three  theories  are  virtually 
identical.  They  arise  from  the  first  antisymetric  branch  motion  of  the  plate 
surface,  viz.  simple  plate  flexure.  The  small  notch  at  17°  in  the  plot  using 
exact  plate  theory  marks  the  initial  .development  of  significant  sound  radi¬ 
ation  from  the  first  symmetric  branch  motion  of  the  plate,  which  is  not  built 
into  the  thin  and  thick  plate  equations. 


At  5kHz,  Figure  5,  the  'hump'  at  ©-€8°  marks  the  development  of  coinci- 


dence  lobe  radiation:  thin  plate  theory  is  now  no  longer  suitable.  At  10kHz, 
Figure  6,  there  is  still  good  agreement  between  the  thick  and  exact  plate 
theories,  except  in  the  region  of  ©-17°  where  the  first  symmetric  branch 
motion  of  the  exact  theory  continues  its  development.  At  20kHz,  Figure  7,  the 
increasing  inadequacy  of  thick  plate  theory  is  evident  at  ©*17°  and  at  all 
angles  greater  than  45°.  At  30kHz,  Figure  8,  thick  plate  theory  is  totally 
inadequate,  except  over  a  small  range  of  angles  near  to  the  coincidence  lobe 
at  e-36». 

At  40kHz,  Figure  9,  the  emergence  of  sound  radiation  due  to  the  second 
antisynmetric  branch  motion  of  the  plate  is  evident  at  ©=9°  in  both  thick 
and  exact  plate  theories.  However,  while  the  thick  plate  theory  is  still 
satisfactory  at  the  coincidence  lobe,  ©-35°,  it  does  not  predict  correctly  the 
level  of  the  lobe  at  ©**9° . 


(<3)  Radiated  Pressure . _ Exact  Theory.  Longitudinal  Force 

Figures  10-11  show  the  variation  with  ©  of  the  far-field  pressure  at 
selected  frequencies  for  the  case  of  longitudinal  point  force  excitation  on 
the  lower  surface,  which  is  in  contact  with  a  vacuum.  The  pressure,  which 
varies  in  the  azimuthal  direction  as  cos($),  is  plotted  at  the  angle  $~o°. 

The  sharp  peaks  at  ©->17°  are  dominant.  When  allowance  is  made  for  the 
inadequate  resolution  of  the  peaks,  the  maximum  levels  are  seen  to  be  close  to 
those  levels  that  would  have  been  obtained  in  the  absence  of  the  plate,  viz. 
20.0  x  log10(Pf/c)  +  120.0,  which  range  from  182. 5dB  at  2kHz  to  208. 5dB  at 
40kHz.  The  levels  of  the  broader  coincidence  lobes,  which  occur  from  lOkKz, 
are  only  3-6dB  less  than  the  corresponding  lobes  due  to  vertical  excitation. 


(e)  Exact  Theory. _ Excitation  on  Upper  Surface 

Comparisons  have  also  been  made  between  power  and  pressure  plots  obtained 
from  point  force  excitation  on  the  upper  and  lower  surfaces,  respectively,  but 
the  plots  obtained  from  the  former  condition  are  not  included  herein.  Power 
levels  show  only  small  differences,  whereas  the  pressure  levels  may  differ 
considerably  in  the  continuum  background  but  not  too  significantly  in  regions 
close  to  the  resonant  peaks. 


(f)  Reflection  and  Transmission.  Exact  Theory 

Figures  12  and  13  show  the  reflection  and  transmission  coefficients, 
respectively,  as  a  function  of  the  angle  of  incidence,  ©.  The  reflection 
coefficient  is  close  to  OdB,  except  at  dips  which  occur  at  the  coincidence 
angles  of  the  symmetric  and  antisymmetric  branch  motions.  At  20  and  40kHz, 
there  are  well  defined  peaks  in  the  transmission  coefficient  at  the  coinci¬ 
dence  angles.  The  reflection  and  transmission  characteristics  of  fluid- loaded 
plates  is  discussed  in  some  detail  in  the  papers  of  Freedman  [6]  and  Fiorito 
[7]. 


(g)  Perspective  and  Contour  Plots.  Exact  Theory 


Figures  14-17  show  perspective  and  contour  plots  of  sound  radiation 


reflection  and  transmission  in  which  the  lose  factor  was  chosen  as  o.'U.  They 
were  obtained  via  the  university  of  Bradford's  SIMPLEPLOT  graphics  package, 
from  a  128x128  array  of  equispaced  decibel  data.  The  perspective  plots  are 
for  a  user  selected  angle  of  view  of  60°  with  the  horizontal;  their  appearance 
is  enhanced  by  setting  a  dynamic  range  of  40dB,  that  is,  values  less  than  the 
maximum  minus  40dB  are  set  to  this  base  level.  The  perspective  plots  enable  a 
quick  qualitative  assessment  of  the  acoustic  properties  of  the  steel  plate  in 
terms  of  the  two  antisynnetric  and  one  synmetric  branch  motions  of  the  plate. 
The  contour  plots  present  the  same  information  mush  less  vividly,  but  are  more 
helpful  if  precise  magnitudes  are  of  concern. 


4.  COMCUJDING  REMARKS 

Numerical  results  obtained  from  exact  linear  theory  have  been  presented 
for  sound  radiation,  reflection  and  transmission  of  a  steel  plate  in  water. 
Sams  of  the  plots  are  new  to  the  literature  and  should  physical  under¬ 
standing  When  studied  in  conjunction  with  some  of  the  papers  referenced 
herein.  The  comparison  of  sound  radiation  predictions  of  thin,  thick  and 
exact  plate  theories  have  illustrated  the  range  of  applicability  of  the 
approximate  plate  theories.  The  predictions  of  sound  radiation  due  to 
longitudinal  excitation  serve  as  reference  plots  for  use  in  the  future  when 
predictions  obtained  from  more  refined  plate  theories  become  available. 


E.J.  CLEMENT  (HSQ) 
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ACOUSTIC  POWER  OF  POINT  EXCITED  PLATE.  B. 
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FIG. 2  POWER  RADIATED  BY  APPROXIMATE  PLATE  THEORIES 


FIG. 3  POWER  RADIATED  BY  EXACT  PLATE  THEORY 
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FIG, 5  COMPARISON  OF  SOUND  RADIATION  AT  5kHz 
VERTICAL  POINT  FORCE  ON  DONER  SURFACE 
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FIG. 6  COMPARISON  OF  SOUND  RADIATION  AT  10kHz 
VERTICAL  POINT  FORCE  ON  LOWHt  SURFACE 
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FIG. 9  COMPARISON  OF  SOUND  RADIATION  AT  *WkHz 
VERTICAL  POINT  FORCE  ON  LOWES  SURFACE 
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FIG.  10  SOUND  RADIATION  OF  EXACT  THEORY  WITH 
LONGITUDINAL  FORCE  ON  LOWES  SURFACE 
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FIG.  11  SOUND  RADIATION  OF  EXACT  THEORY  WITH 
LONGITUDINAL  FORCE  ON  LOWER  SURFACE 


TRANSMISSION  COEFFICIENT  OF  STEEL  PLATE.  EXACT  THEORY. 


VfflTICAL  EXCITATION.  CQKTOUBS  IN  MULTIPLES  OF  10. 
MAXIMUM  LEVEL  -  205.(5  dB.  MINIMUM  LEVEL  -  113.4  ctB. 
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